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Abstract 

We consider the monopole-antimonopole static potential in the confining phase of the 
Abelian Higgs model and in particular the corrections to the Coulomb-like potential at 
small distances r. By minimizing numerically the classical energy functional we observe a 
linear in r, stringy correction even at distances much smaller than the apparent physical 
scales. We argue that this term is a manifestation of the condition that the monopoles 
are connected by a mathematically thin line along which the scalar field vanishes. These 
short strings modify the operator product expansion as well. Implications for QCD are 
discussed. 



"String" or vortex-like configurations play an important role in a number of areas of physics 
and in theoretical speculations. The best known perhaps is the flux tube of superconductiv- 
ity ("Abrikosov string") and its relativistic version in particle physics, in the Abelian Higgs 
model ("Nielson-Olesen string") 0. Such configurations are also believed to play an impor- 
tant role in QCD, and the confining force between two quarks is often described as due to the 
stretching of the string 0. 

Such physical strings have a definite thickness, reflecting the balance of the various forces 
going into their composition, and this introduces a certain length scale into the problem. For 
example for the forces between two static sources, like the quarks of QCD or between hypo- 
thetical magnetic monopoles in a superconductor, the string gives a linear potential energy at 
large distances, reflecting a certain energy per unit length of the string. 

Naturally at small distances, that is distances small compared to the thickness of the string, 
where the picture of an energy per unit length does not seem to apply, we would expect the term 
linear in the separation to disappear. In this paper we would like to point that, surprisingly, 
even at small distances, that is less than the thickness of the string, there remains a stringy 
term in the energy. This term may be interpreted as physical manifestation of the mathematical 
Dirac string which accompanies such objects as magnetic monopoles. 

In the bulk of the paper we will consider the Abelian Higgs model (AHM) while implications 
for QCD are summarized in the concluding remarks. The AHM describes a gauge field 
interacting with a charged scalar field $ as well as self-interactions of the scalar field, and the 
corresponding action is: 

S = jd*x {±1* + \\{d - zA)$| 2 + iA(|$| 2 - V 2 ) 2 ) (1) 

where e is the electric charge, A, rj are constants and F^ v is the electromagnetic field-strength 
tensor, F^ u = d^A v — d u A^. The scalar field condenses in the vacuum, (<£>) = rj and the physical 
vector and scalar particles are massive, my = e 2 r] 2 ,m 2 H = 2Xrj 2 . 

The model is famous to provide with a relativistic analog of the Landau- Ginzburg theory of 
superconductivity. In particular, if one introduces a monopole-antimonopole pair as an external 
probe its static potential V(r) grows linearly with the distance r at large r: 

Urrir^ooVir) = a^r (2) 

The growth of the potential (2) is well understood in terms of the Abrikosov-Nielsen-Olesen 
(ANO) strings Q. These strings are solutions to the classical equations of motion corresponding 
to the action (1) and carry (quantized) magnetic flux equal to 27in/e. Because of this, strings 
may end up with monopole-antimonopole pairs. The salient features of the solution is that the 
Higgs field $ disappears on the axis of the string and reaches its vacuum value at transverse 
distances of order l/m# while the magnetic field extends to distances of order 1/my. In the 
limit of large distances, r 3> m v 1 1 m'J^ the ANO string can be considered as thin and the 
constant in Eq. (0) is then its tension. 
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We will study the static energy of a monopole-antimonopole pair at short distances, r <C 
rriy ,m]j. The potential V(r) is then Coulomb like because of exchange of the photon. As is 
mentioned above, the ANO strings are not relevant at such distances, and one would expect 
that the power corrections to the Coulomb potential are simply due to a non-vanishing vector 
particle mass. That is, up to a constant which can be included into definition of the heavy 
masses: 



Note the change of the sign in front of the linear term as compared to the stringy potential (§). 

The central point of the present letter is that these naive expectations are not true because, 
even at scales when the ANO string is irrelevant, there still exists an infinitely thin line inherent 
to the problem. Namely, there is a line connecting the monopoles along which $ = 0. Its 
existence follows from the observation || that the world sheets with $ = are either closed 
or have monopole trajectories as their boundaries. In other words, it is a topological condition 
that magnetic charges are connected by a line $ = 0, no matter how small the distance r 
between the charges is. 

Because this topological condition is so important for further results, it is worth noting 
that it can be rederived in the language of the Dirac string || as well. Namely, the infinitely 
thin line discussed above is nothing else but the Dirac string connecting the monopoles. The 
possibility of its dynamical manifestations arises from the fact that the Dirac string cannot 
coexist with $ ^ and $ vanishes along the string. Indeed, the self-energy of the Dirac string, 
is normalized to be zero in the perturbative vacuum. To justify this one can invoke duality and 
ask for equality of self-energies of electric and magnetic charges. Since the electric charge has 
no string attached the requirement would imply vanishing energy for the Dirac string. However, 
if the Dirac string would be embedded into a vacuum with ($) ^ then its energy would jump 
to infinity since there is the term 1/2|$| 2 A^ in the action and A 2 ^ — > oo for a Dirac string. 
Hence, $ = along the string and it is just the condition found in |3J. In other words, Dirac 
strings always rest on the perturbative vacuum which is defined as a vacuum state obeying the 
duality principle. Therefore, even in the limit r — > there is a deep well in the profile of the 
Higgs field This might cost energy which is linear with r even at small r. 

By solving numerically the classical equations of motion with a boundary condition $ = 
along the line connecting the monopoles we do find a linear stringy piece, i.e. with a positive 
coefficient in front of r, in the potential even in the limit r — > 0. Another manifestation of the 
fact that we are dealing with a kind of elementary string is the breaking of the operator product 
expansion. Indeed the standard operator product expansion (OPE) assumes that short distance 
expansions can be derived in terms of exchange of elementary particles. On the other hand, 
the topolgical string $ = cannot be constructed in terms of particle exchanges but should 
be added independently. Since the effective string tension for small size strings is proportional 
to < $ > 2 , the standard OPE breaks down at the level of r 2 < $ > 2 corrections to pure 
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perturbative results. 

To evaluate the static energy E(r), we first need to introduce the monopole trajectories 
explicitly into the path integral: 

E(r)=m^(~ > jhiH(C) t (3) 

H{C) = \\ DAD$[D$}H(A, E c )e~ S{A ^ , (4) 
here the 't Hooft loop is defined as follows: 

H(A, £ c ) = exp J d * x _ + 2vr [Z c ]y] } (5) 

where S c is an arbitrary surface having the contour C as a boundary, 5T, C = C . A particular 
form of the contour C is determined by the monopoles trajectory and for a static monopole- 
antimonopole pair located at distance r the contour C is a rectangular loop T x r with T 3> r. 
In Eq (H) [...] denotes the duality operation so that for any tensor Tp V of the second rank 
[T/u] = \^^u\p T\ p . It can be shown that the expectation value of the 't Hooft loop does 
not depend on the particular choice of the surface £ . Moreover it can be shown JH] that 
in the string representation of the Abelian Higgs model we have to sum over all surfaces oq 
spanned on the loop C and also over all closed surfaces a (representing virtual glueballs). The 
classical solution corresponds to the minimal area surface oc- Here we consider the classical 
approximation, hence the problem of finding Eir) is equivalent to solving classical equations 
of motion with boundary condition $ = along the straight line connecting the monopoles. 

Let us note that the numerical results for the potential V(r) can be found in a number 
of papers |6], 0] and in this respect our approach is not new. However, there are no measure- 
ments dedicated specifically to small corrections to the Coulombic potential at r — > 0. Indeed 
the standard question addressed is how the Coulomb-like behaviour of the potential at short 
distances is transfromed into the linear potential at large distances (see, e.g., Ref. || and 
references therein). It is only very recently that it was recongnized that the power correction 
to the potential at short distances in QCD could signify a new phyisics 0. Hence, there are no 
measurements which would provide error bars of the slope of the potential at small distances. 
Moreover, one could argue that at small r the strong magnetic field "pushes out" the Higgs 
field in any case and, therefore, the potential energy is insensitive to the condition $ = along 
a line connecting the magnetic charges. Since in answering this kind of questions we rely on 
the numerical results, we need dedicated measurements. 

We will consider the unitary gauge, Jm$ = 0. Then the most general ansatz for the fields 
consistent with the symmetries of the problem is: 

$ = r]f (p, z) A a = £ ab x b A(p, z) A = A 3 = 

(6) 

P = [x a x a } 1/2 z = x 3 x a = x a /p a = 1, 2 . 
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Let us introduce also a new variable k = v^2A/e = m H /m v and measure all dimensional 
quantities in terms of m v = erj. Then the energy functional is: 



E{r) = £ / dz J pdp{[\d p (pA) + £] 2 + [d z Af+ 



+ [d p /] 2 + [d 2 /] 2 + PA* + \k\P - l) 2 } 



1 

Z = -5(p)- J d£6(z-Q 



(7) 



(8) 



In the limit r — > the Coulombic contribution becomes singular. The easiest way to separate 
the singular piece is to change the variables A = A^ + a, where Ad is the solution in the absence 
of the Higgs field: 



1 

2p 



r + 



z± = z ± r/2 r ± 



2 i 2 



1/2 



(9) 



Upon this change of variables the energy functional takes the form: 
E(r) = E self -n/r + E(r), 

E(r) = ^ J dz J pdp{[%(pa)f + \d z af+ 

-oo k 1 

+ [d P f] 2 + [d z ff + P(a + A d f + \k\P - l) 2 } 



(10) 



We investigated numerically E{r) in the region 0.1 < r < 5.0 (details of the procedure 
will be given elsewhere). Fig. 1 represents the results of our computations for various k values 
and clearly demonstrates that there is a linear piece in the potential even in the limit r C 1. 
While obtaining the values of the slope <jq at small distances we separated the Yukawa-potential 
contribution by fitting the energy as follows: 



E fit (r) = C 



l-e 



1 )+(cj + -C )r = a r + O(r 2 ) 



(11) 



The slope a depends smoothly on the value of k, see Fig. 2. The linear piece in the potential 
at small distances reflects the boundary condition that $ = along the straight line connecting 
the monopoles. To illustrate this point we show in Fig. 3 the function (1 — f(p, z)) which makes 
the impact of the stringy boundary condition visible. It is noteworthy that in the Bogomolny 
limit (k—1) the slope of the linear potential within the error bars is the same for r — > oo and 
r -> 0. 
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As it is mentioned above, the existence of short strings is manifested also through breaking 
of the standard operator expansion. Indeed, above we found the potential in the classical 
approximation. In this approximation the potential is usually directly related to the propagator 
D^ u (q 2 ) in the momentum space, 

V(r) = jtfVe tq - r A)o(q 2 ) (12) 

Moreover, as far as q 2 is in Euclidean region and much larger than the mass parameters, the 
propagator D^^q 2 ) can be evaluated by using the OPE. Restriction to the classical approxi- 
mation implies that loop contributions are not included. However, vacuum fields which are soft 
on the scale of q 2 can be consistently accounted for in this way (for a review of see, e.g., ||10||). 
This standard logic can be illustrated by an example of the photon propagator connecting two 
electric currents. Modulus longitudinal terms, we have: 




\ ,,_ gz qz qz qz qz j qz _ m ^ 



Thus, one uses first the general OPE assuming \q 2 \ ^> e 2 $ 2 then substitutes the vacuum 
expectation of the Higgs field $ and upon summation of the whole series of the power corrections 
reproduces the propagator of a massive particle. The latter can also be obtained by solving 
directly the classical equations of motion. 

This approach fails, however, if there are both magnetic and electric charges present. In this 
case, one can choose the Zwanziger formalism [llj and work out an expression for propagation 



of a photon coupled to magnetic currents following literally the same steps as in (^3|), (see e.g. 
31). In the gauge n^D^ 




D, v (q,n) = ^^^-^(^ + ^) + ^ + ^(^n 2 -n M n,)J. (14) 

Here the vector n M is directed along the Dirac strings attached to the magnetic charges and 
there are general arguments that there should be no dependence of physical effects on |T 



On the other hand, if the potential energy is given by the Fourier transform of (II) then its 
dependence on n M is explicit, see, e.g., |L3j, [L4j and we addressed this problem in Ref. IJ. 

Note that Eq. flU] ) immediately implies that the standard OPE does not work any longer 
on the level of q~ 2 corrections. Indeed, choosing q 2 large and negative does not guarantee now 
that the rriy correction is small since the factor (qn) 2 in the denominator may become zero. 
Of course, appearance of the poles in (qn) variables in longitudinal pieces is not dangerous 
since they drop due to the current conservation. However, the term proportional to rriy in ( |l4l) 
cannot be disregarded and contribute, in particular, to ((T2j) . 

The reason for the breaking of the standard OPE is that even at short distances the dy- 
namics of the short strings should be accounted for explicitly. In particular, in the classical 
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approximation the string lies along the straight line connecting the magnetic charges and af- 
fects the solution through the corresponding boundary condition, see above. More generally, 
the OPE allows to account for effect of vacuum fields, in our case for ($) ^ 0. The OPE is 
valid therefore as far as the probe particles do not change the vacuum fields drastically and the 
unperturbed vacuum fields are a reasonable zero-order approximation. In our case, however, 
the Higgs field is brought down to zero along the string and this is a nonperturbative effect. 
Thus, the stringy piece in the potential V(r) at r — > is a nonperturbative correction which is 
associated with short distances and emerges already on the classical level. 

A few remarks on the implications of the results obtained to QCD are now in order. To 
begin with, there exist detailed numerical simulations on the lattice which confirm the dual- 
superconductor picture of the confinement (for review and references see []15|). In particular, 
the potential V(r) at large r is well described by the model ||. Moreover, if the simulations are 
performed in the U(l) projection of QCD, condensation of a scalar field <3? ma g with magnetic 
charge is confirmed and, moreover, the structure of the observed string which determines the 
QQ potential at large distances is well described by the classical Landau-Ginzburg equations 
||16|| . What is especially important for us, is that the definitions of the magnetic charges in 
the Abelian projection of gluodynamics Q are in fact local and, therefore, the results obtianed 
within the AHM with a local field $ can imitate gluodynamics. 

Note that, naively, existence of a U(l) gauge invariant operator |$ ma9 | 2 of dimension d = 2 
in the abelian projection of QCD would imply infrared sensitive corrections of order ($ mas ) 2 /<? 2 
which are calculable via the OPE. On the other hand, it is well known (for references see, 
e.g., [[[(J) that such corrections are not present in QCD. The paradox is resolved through the 
observation that the U(l) projection of QCD is similar to the Higgs model and, therefore, the 
OPE breaks down in this projection. 

Moreover, we have learned that in the Higgs model there emerges in fact a non-perturbative 
short distance correction of order q~ 2 manifested through the slope do 7^ 0. Although the Dirac 
strings are specific for the abelian projection, the results for the slope <Jq which is a physical 
quantity should be true for any gauge fixing. In other projections, the linear potential at small 
distances arises if there are short- distance correlations in non-perturbative fields ||. There 
are attempts to extract phenomenological consequences from hypothetical existence of short 
strings in QCD f20 |. On the other hand, the most common picture of non-correlated finite size 



non-perturbative fluctuations results in cr = flT7|. 

It is amusing therefore that the lattice simulation [IS] do not show any change in the slope 
of the QQ potential as the distances change from the largest to the smallest ones available. In 
the notations introduced above, 

o-Qc ~ <Tq. (15) 
Moreover, it is known from phenomenological analisis and from the calculations on the lattice 



13| , 18| , |16f that the realistic QCD corresponds to the f=a 1 where k = vein fray. It is 
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remarkable that, as is mentioned above, the AHM in the classical approximation also results 
in the relation (|15[) for /c»l. Thus we see that existing data |Tj| on the behavior of the QQ 
potential at small distances agree with the classical approximation to AHM. Also, our results 
support indirectly phenomenological attempts to account for the novel 1/q 2 corrections [^0 



To summarize, we have demonstrated that the potential of a monopole-antimonopole pair 
at distances much smaller than the inverse masses rriy l H does contain a linear piece cr r with a 
positive cr - This linear piece is a dynamical manifestation of the topological condition that the 
scalar field <3> vanishes along a line connecting the magnetic charges. These short strings are 
responsible also for breaking of the standard OPE on the level of 1/q 2 corrections. Note that 
usually the Abelian Higgs model plays the role of the effective infrared model of gluodynamics 
(see reviews [15]). It is amusing that the behaviour of the QQ potential at small distances 
obtained via the lattice simulations agree with the dual-superconductor model with 
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Figure 2: The linear slop a of the energy ^ • ^- in the limit r — > as a function of k (see Qll 
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Figure 3: The function (1 — f(p, z)) in the p — z plane for ran = f^v — m, r — 0.2/m. The 
line at which / = is clearly seen. 
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